Predicting B° -> p°j and B° -> p°7 using holographic AdS/QCD 
Distribution Amplitudes for the p meson 



M. Ahmady 

Department of Physics, Mount Allison University, Sackville, N-B. E46 1E6, Canada* 

R. Sandapen 
Departement de Physique et d'Astronomie, 
Universite de Moncton, Moncton, N-B. E1A 3E9, Canada & 
Department of Physics, Mount Allison University, Sackville, N-B. E46 1E6, Canada t 

We derive holographic AdS/QCD Distribution Amplitudes for the transversely 
polarised p meson and we use them to predict the branching ratio for the decays 
B° — > p°7 and B° — > p°j beyond leading power accuracy in the heavy quark limit. 
For B° — > p°7, our predictions agree with those generated using Sum Rules (SR) 
Distribution Amplitudes and with the data from the BaBar and Belle collaborations. 
In computing the weak annihilation amplitude which is power-supressed in B° — > p°7 
but is the leading contribution in B° — > p°7, we find that, in its present form, the 
AdS/QCD DA avoids the end-point divergences encountered with the SR DA . 

INTRODUCTION 

Flavor Changing Neutral Currents (FCNC) are excellent probes to the Standard Model 
and beyond. In particular, the b — > (s, d)j transitions are most important for the extraction 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements as well as for the search of 
New Physics (NP) signals. The experimental measurements for the b — > d^y transistion are 
currently available for exclusive radiative B decay to a p meson, i.e. B — > pj. For a recent 
review of radiative B decays, we refer to [1]. 

The theory of exclusive decays is complicated by their sensitivity to non-perturbative 
physics. The standard[40] theoretical framework is QCD factorization (QCDF)[2, 3] which 
states that, to leading power accuracy in the heavy quark limit, the decay amplitude factor- 
izes into perturbatively computable kernels and non-perturbative objects namely the B — > V 
transition form factor, the meson couplings and the leading twist Distribution Amplitudes 
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(DAs) of the mesons. Traditionally the DAs for the vector meson are obtained from QCD 
Sum Rules[4-6]. The numerical values of the transition form factor and the tensor coupling 
of the vector meson are obtained from light-cone Sum Rules or lattice QCD. The predic- 
tive power of QCDF is limited by the uncertainties associated with these non-perturbative 
quantities and also by power corrections to the leading amplitude [7]. The computation of 
the power corrections is often problematic due to the appearance of end-point divergences 
in convolution integrals that contribute to the decay amplitude [7, 8]. 

Our goal in this paper is to use new holographic AdS/QCD DAs for the p meson to 
compute the branching ratios for two exclusive radiative decays, namely B° — > p°j and B° — > 
p°7, beyond leading power accuracy. We derive the AdS/QCD DAs using a holographic 
AdS/QCD light-front wavefunction for the p meson [9] which was recently shown to generate 
predictions for the cross-sections in diffractive p meson production that are in agreement 
with the data collected at the HERA electron- proton collider [10]. Reference [10] also shows 
that the second moment of the AdS/QCD twist- 2 DA of the longitudinally polarised p 
meson is in agreement with both Sum Rules and lattice predictions. Here we shall extend 
the comparison between AdS/QCD and Sum Rules for the twist- 2 and twist-3 DAs of the 
transversely polarized p meson since both of these non-perturbative quantities are required 
to compute the decay amplitudes for B° — > p°7 and B° — > p°7 beyond leading power 
accuracy. 

Theoretical predictions for the branching ratio of B° — > p°7, using the leading twist-2 SR 
DA for the p meson, can be found in references [2, 11]. In both references, the leading power 
correction (0(A QCD /m5)) due to annihilation is taken into account within QCDF and the 
subleading (0(A QCD /mb) 2 ) annihilation contributions are neglected. On the other hand, in 
reference [11], two other classes of power corrections due to long-distance photon emission 
and soft gluon emission are also taken into account. Here we shall investigate the numerical 
importance of three additional subleading annihilation contributions to B° — > p°7, two of 
which turn out to be sensitive to the higher twist-3 DA of the p meson. In fact, the rare 
decay B° — > p°7 proceeds mainly via these four annihilation processes and a theoretical 
prediction for its branching ratio using the SR twist-3 DA is available in reference [12]. 
However this prediction suffers from a large uncertainty because of end-point divergences 
encountered when computing those annihilation contributions sensitive to the twist-3 DA 
for the p meson. We shall update this prediction using the twist-3 AdS/QCD DA which, as 
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we shall see, avoids the end-point divergence problem. 

On the experimental side, the branching ratio for B° — > p°j has been measured with 
increasing precision by the BaBar and Belle collaborations [13]. On the other hand, the rare 
decay B° — » p°7 has not been measured experimentally but it is an interesting process to 
investigate at LHCb because of its sensitivity to NP especially those which allow FCNC at 
tree level [12]. 

We start by an outline of the computation of the amplitudes for both decays following 
references [2, 12], devoting attention to the quantities which are dependent on the DAs of 
the p meson. 

DECAY AMPLITUDES 

The effective weak Hamiltonian for the underlying b — > dry transition is given by 

U ^ = % E v pM c iQi + c *q p 2 + E w*i w 

p=u,c i=3 

where q = d for B° — > p 0/ ~f and q = s for B° — > p 0/ ~f. Q\ and Q\ are the current-current 
operators, Qi=s..e are the QCD penguin operators and Q7 and Q$ are the electromagnetic 
and chromomagnetic operators. The coeffecients Cj are the perturbatively known Wilson 
coeffecients and are the CKM matrix elements. In this paper, we shall use the NLO 
Wilson coeffecients given in [11, 14] [41] and we use the numerical values of CKM matrix 
elements given in reference [15]. 

We start with the amplitude for the decay B° — > p 0/ ~f. At leading power accuracy in the 
heavy quark limit and to all orders in the strong coupling a s , the matrix elements of these 
operators factorize as [2]: 

(p(P,e T Mq,e)U\B) = [F B ^T! + f d( dz <& B (0?f '(C, • * (2) 

Jo 

i.e. into perturbatively computable hard-scattering kernels T/ and T t n and three non- 
perturbative quantities namely the transition form factor F B ^ P , the leading twist DA of 
the B meson, $s(C), and the twist-2 DA of p meson, (fif(z). In equation (2), P and ex are 
the 4-momentum and polarization vector of the p meson while q and e are the 4-momentum 
and polarization vector of the photon. The form factor F B ^ P is obtained from light-cone 
QCD Sum Rules [11]. The second term in equation (2) describes mechanisms involving 
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the spectator quark, hence its dependence on the DAs of the mesons. In what follows, we 
shall only need the first inverse moment of $b(C) which is parametrized as M B /\ B where 
\ B = 0(A QCD ) [2]. 

At zeroth order in a s , the leading power amplitude is given by 

^ p g o 7 = %£ V; d V pb C 7 (Q 7 ) (3) 

p=u,c 



where (Q 7 ) is the matrix element of the operator Q 7 , i.e. 

i b {p)F 

2V27T 



(p(P, eh(q, e)\Q 7 \B°) = [ we ^pV + i((c . eT )( P ■ g) - ( e ■ P)(g • e T ))] 



(4) 

where e = \/47rQ! cm and m b (p) is the running mass of the b quark evaluated at the hard scale 
p = m b . 

At order a s , the leading power amplitude becomes [2] 

Gt 



^ a :; S o 7 = J E ^«?<0r> (5) 



with 



a? = C' 7 + ^^[C 1 ( /U )G 1 ( Sp ) + C 8 ( / i)G 8 ] 

+ ^[CiWli( Spi ^) + C 8 W// 8 W] . (6) 

47T 

In equation (6), the strong coupling is evaluated at two different scales: \i = m b and a 
hadronic scale /j,h = \J A-qcd/J'- The functions G\ and Hi depend on s p = (m p /m b ) 2 where 
m p is the quark mass in loops contributing at next-to-leading order accuracy in a s . The 
hard scattering functions Gi(s p ) and Gg are given explicitly in reference [2]. Here, we focus 
on the functions Hs(/J>h) and Hi(s p , jj, b ) which depend on the twist-2 DA of the p meson. 
The function Hi(s p , p) is given by [2] 

where f B is the decay constant of the B meson which is obtained from lattice QCD [16, 17]. 
i] w2 (s p , p) is a convolution of the twist-2 DA with a hard scattering kernel, i.e. 



IT\s p ,p)= [ dzh(s p , 
Jo 



z)<f>j(z,p) (8) 
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where the hard scattering kernel is given by 

4Sr 



h(s p ,z) = 



z 2 



+ L 2 



(9) 



with L 2 being the dilogarithmic function and z — 1 — z. The function H 8 (ii) is given by 

;) 



HM ~ I 37V C F^PM1 J ) h (10) 



where /^(aO i s the first inverse moment of the twist- 2 DA, i.e. 



-(z,n) 



ITUi) = / dz . (11) 
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We note that if m& ^> m p , then 

/r 2 ( Sp ,/i)^/r 2 (o,/i) = -2/r(^) (12) 

so that both Hi and H 8 become simply proportional to the first inverse moment of the 
twist-2 DA of the p meson. In practice, this approximation is not justified for a charm loop, 
i.e. when p = c and we do not make it here. In what follows, we shall take m u ^ = 0.14 GeV 
[10], m c = 1.3 GeV and m b = 4.2 GeV. 

In B° — > p°7, all annihilation topologies are suppressed by at least one power of A QCD /m b 
[2]. The leading annihilation contribution is given by [2] 

A^LuoSB P°7) = ^§Ok (c 2 + b d (Q 7 ) (13) 

with 

_ 2Kf B f p M p 

bd ~ FB-+PM*\ B ■ (14) 
This leading contribution can be taken into account by adding an extra term to the coeffe- 
cient ay in the leading power amplitude given by equation (5): 

a u^ a u + bd ( C2+ _L C \ (15) 



where bd is given by equation (14). The leading annihilation contribution corresponds to the 
annihilation diagram in which the photon is radiated off the spectator quark of the B meson, 
i.e. the third diagram in figure 1. Here we wish to investigate the numerical importance 
of the three other subleading annihilation contributions shown in figure 1. In fact, the four 
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annihilation diagrams of figure 1 are the dominant contributions to the decay B° — > p°7 
[12]. The total annihilation amplitude is given by [12] 

eG Fl 



^annihilation _ _____ y*^ ^ V tb f B(3) f pM p {A A ( s ) + A^) + A^ + A\( g ) ) 



where to zeroth order in a. 



(16) 



(17) 



with 



and 



while 



rtw3 



m = r 

^0 



#p w (z,/z) 



dz 



V(s) P J 

Mv), 



1 dz zgp {V, (z,fi) 



o 



zMl, + zzM„ 2 - ml 



f 



A, 



+i 



2(e.P)(_.e T )-M^ (s) 1 



M flw 



(e • e T ) 



and 



m 6 



2(e.P)(q.e T )-M' B(s) M + -^l( c . eT ) 



(18) 



(19) 



(_^^ ) ( W4^ Q / + <[(c • er)(P • ?) - (c • • er)] (20) 



{ W e " e T*V - • er)(P ■ q) - (e ■ P)(q ■ e T )\ (21) 



In the above equations, C12 and C 34 are the combinations of the Wilson coefficients [42] [12]: 
1 



C 



12 



C34 — 



v/2 
1 

71 



2 / ^\ / ^ ^ 

3 V 3 ) V tb V t * d{s) ^ V 3 3 

~ CA K&Kx*_( s ) 

3 v 2 3 y \ ;/,\ ; 



td(s) 



(22) 
(23) 



and E~f 8 \ is the energy of the photon in the Bu\ meson rest frame, i.e. 



-& 7 ( s ) = — z — I 1 



ME 



(24) 
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FIG. 1: Annihilation processes contributing to the decay B° — > p°^f and B° — > p°j. The diagrams 
(1),(2) and (4) are power-suppressed compared to the diagram (3). The contributions (1) and (2) 
are sensitive to the twist-3 DA of the p meson. 

The quantity \b s is analogous to A#, i.e. it parametrizes the first inverse moment of the 
B° meson DA. As expected the amplitudes A^ s y corresponding to annihilation topologies 
in which the photon is radiated off the light quark or antiquark of the p, are sensitive to 
the twist-3 DA of the p meson. Note that the annihilation amplitude A% evaluated to 
leading power in the heavy quark limit, coincides with equation (13) which is the leading 
annihilation contribution given in reference [2]. 
The total decay amplitude for B° — >■ p°j is then 

A(B° -> p° 7 ) = ^" p g o 7 + A^ nmatioa (25) 

where *4.^ c a ^ ^ and ^4. AnniliUatlon are gj ven equation (5) and (16) respectively. On the other 
hand, the total decay amplitude for B° — > p°7 is given by 

A(B° p°7) = A? nihilation (26) 

where ^Annihilation j g gj ven j^y equation (16). 

In order to compute the decay amplitudes given by equations (25) and (26), we must 
specify the twist-2 DA 0^ in equations (8) and (11) as well as the twist-3 DA g 1 ^ appearing 
in equations (18) and (19). We also need to specify the numerical value of the tensor coupling 
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fp which appears in equations (7) and (10). We shall do this in the next section using a 
holographic AdS/QCD light-front wavefunction for the p meson. The numerical values of 
the decay constants /b (s) , the parameters A( s ) and the form factor F B ^ p are shown in table 
I. 



Parameter 


Numerical value 




0.51 ±0.12 (0.60 ±0.20) GeV 




190.6 ± 47 (227.6 ± 5.0) MeV 


pB^p 


0.27 ±0.04 



TABLE I: The non-perturbative input parameters obtained from QCD Sum Rules or lattice QCD 
[11, 16-19]. We shall take the central values of these parameters to compute our predictions for 
the branching ratios. 



HOLOGRAPHIC ADS/QCD DAS AND COUPLINGS OF THE p MESON 

Distribution Amplitudes parameterize the operator product expansion of vacuum-to- 
meson transition matrix elements of quark- ant iquark non-local gauge invariant operators 
at light-like separations. At equal light-front time x + = and in the light-front gauge 
A + = 0, we have [4, 5] 

<0|9(0)y*gOO|A>(P,A)> = f p M p ^P» j\ue- iuP+ *-<t>l{u,p) 

+ f p M p (< - ^ J^) f du e-^V'W) , (27) 

(0\q(0)kf,'f\q{x-MP,\)) = 2ft{£P"-e v x P'>) f du e~ iuP+x ~ ^ («, v) (28) 
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and 



(0\q(0)ri 5 Q(x-)\p(P,\)) = --^ pa e\PWfpMpj\u ^^(u, p) (29) 

for the vector, tensor and axial-vector current respectively. The polarization vectors e\ are 
chosen as 

/ P + M \ 1 
eL= {w~P+' 0± ) e r(±) = -^(0,0,l,±i) (30) 
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where P + is the "plus" component of the 4-momentum of the meson given by 

/ M 2 \ 

P^[P + ^,0 ± j . (31) 



All four DAs satisfy the normalization condition 

/ dz V (z,p) = l (32) 
Jo 

where <p = {4>p' ± , gp^'^} so that for a vanishing light-front distance x~ = 0, the definitions 
of the vector coupling f p and tensor coupling are recovered, i.e. 

(0|g(0) 7 ^(0)|p(P,A)) = / p M^ (33) 

and 

(0|g(0)[ 7 ^7lg(0)|p(P, A)) = 2#(e£P" - e'P") . (34) 

The vector coupling f p is accessible experimentally via the leptonic decay width of the p 
meson [11] 

where T e + e - = 7.04 ± 0.06 keV [15]. On the other hand, the tensor coupling is not 
measured experimentally but is predicted theoretically by QCD Sum Rules and lattice QCD. 
It follows from equations (27), (28) and (29) that the twist- 2 DAs are given by 

f p <f>l(z,p) = J dx~ e- p+ -(0|g(0) 7 + g(x-)|p(P,L)) (36) 



and 
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fpWw) = 2 J dx ~ e X ~mOW n± yl,J + ]q(x-)\p(P,T(±))) (37) 
while the twist-3 DAs are given by 

f P 9p {v) (z,p) = ^Jte~ e- P+ ^(0|g(0)e^ (±) .7g(x-)|p(P,T(±))) (38) 



and 

.(a) 



c\a w 2P+ f 

U-^-M = ^ J dx ~ e l2P+ ^(0|g(0)e^ (±) .77 5 g(^)|p(PT(±))) . (39) 

To relate the DAs to the light-front wavefunctions of the p meson, we use the relation 
[20] 

/at r\k\<n j2u 

dx-e*-* p+ (0|g(0)rg(x-)|p(P,A)) = ^ £ J j-y 2 S$(z, k)0 A (z, k) 

h,h 

, v- h ((l-z)P+,-k) T u h (zP+,k) 
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*)\(z,() = N\-^=^/ z{\ - z) exp f —J exp 



rrif 

(41) 



where we have identified the renormalization scale // as a cut-off on the transverse momentum 
of the quark [20] and A (z,k) is the meson light-front wavefunction in momentum space. 
A two-dimensional Fourier transform of <f>\(z, k) gives the light-front wavefunction, (f>\(z, r), 
in configuration space. The light-front wavefunctions can be modelled [21-23] or extracted 
from data [20]. Here we use the AdS/QCD holographic wavefunction predicted in [24, 25] 
and which can be written as [26] 

i 

2k 2 z{1-z) 

where ( = \/ z(l — z)r is the transverse distance between the quark and antiquark at equal 
light-front time [43] and is the variable that maps onto the fifth dimension of AdS space 
[9, 27, 28]. The AdS/QCD wavefunction given by equation (41) is obtained using the 
soft-wall model [29] to simulate confinement and in that case the parameter k = M p /\^2 
where M p is the mass of the p meson. This AdS/QCD wavefunction has recently been 
used within the dipole model to generate parameter-free [44] predictions for diffractive p 
meson electroproduction that are in agreement with the HERA data [10]. As discussed in 
reference [10], the normalization A/a of the AdS/QCD wavefunction is allowed to depend on 
the polarisation of the meson A = L,T. 

Going back to equation (40), the spinor wavefunctions S^(z, k) are given by [20] 



^k) = 



M, 



p 



m 2 + k 2 

*(1 - z)M 2 



k,-h (42) 



and 



S P h l (± \zM) = ^^y{[(l - z)5 h ^- h± - z5 h± ^\ke ±ie « =F m f 5 h± - h± } (43) 

while T stands for 7 + , [e^^yj, 7 + ], e* T( ^ ± y^ or The matrix element in curly 

brackets of equation (40) can then be evaluated explicitly for each case using the light-front 
spinors of reference [30]: 

" 7 + ^ = 2P+V^ (44) 



^==[4( ± )-7,7 + ]^| = ^V2P+5 h±M , (45) 

^7=4 (±) -7^| = -7^-Al(l - *)<W ± - zS^ke^* =F m f 6 h± - h± } (46) 
V(l - z) y ' \/ z z{l - z) 
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and 



^=4( ±) -77 5 ^| = -^—AT[zS h±rhT + (1 - z)5 hTrh± ]ke^ + (1 - 2z)m f S h± - h± } 

V(l - z) V z z ( l - z ) 



(47) 



We are then able to deduce that 



^,,) = ^/drMW^|, (49) 

9 " < " ,( -" ) = 2^ / dr K J ^ r) W " (z2 + (1 " ^ )2)V - 2 ] (50) 



and 



-(«) 



d, ^") = S/ drM0,r)(1 - !!z)[m '-^^$- (51) 

Equations (48) and (50) were derived in reference [20] where the light-front wavefunctions 
4>\(r,z) were extracted from data. Equations (49) and (51) are new results which show 
how the twist-2 and twist-3 DAs of the transversely polarised p meson are related to its 
light-front wavefunction. 

We are also able to express the vector and tensor couplings f p and in terms of the 
light-front wavefunctions. From the definitions (33) and (34), it follows that 

(0\q(0y L -iq(0)\p(P,L)=f p M p (52) 

and 

<0|g(0)[4 (±) - 7 ,7 + ]g(0)|p(P,T)) = 2f^P + . (53) 

After expanding the left-hand-sides of equations (52) and (53), we obtain the decay width 
constraint [22] 



(54) 

r=0 



and 



n Jo J z(L- z) 

respectively. Note that equations (54) and (55) can also be obtained by inserting equations 

(48) and (49) into the normalization conditions on the twist-2 DAs, i.e. into 

dz <fA(z, oo) = 1 (56) 



/ 

Jo 
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and 

&z<ft{z,n) = l (57) 

'o 

respectively. 
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COMPARISON TO DAS AND COUPLINGS FROM SUM RULES 



Inserting equation (41) in equations (54) and (55), we can compute the AdS/QCD pre- 
dictions for the vector and tensor couplings of the p meson. We compare our predictions 
to experiment, Sum Rules and the lattice in table II. As can be seen, there is reasonable 
agreement between the AdS/QCD prediction for the vector coupling f p and experiment. 
We note that our prediction for fp(fj) hardly depends on /i for /i > 1 GeV. Viewed as a 
prediction at some low scale \i ~ 1 GeV, the agreement with Sum Rules and the lattice is 
reasonable. We note that the AdS/QCD prediction for the ratio of couplings is sensitive to 
the quark mass. For instance, using a current quark mass would yield a ratio far lower than 
the Sum Rules and lattice predictions. We therefore use here a constituent quark mass of 
0.14 GeV which is also the value used in [31-33]. 

Couplings of the p meson 



Reference 


Approach 


Scale p 


fp [MeV] 


fpiv) [MeV] 


frM/fp 


[15] 


Experiment 




220 ±2 






This paper 


AdS/QCD 


~ 1 GeV 


214 


135 


0.63 


[11] 


Sum Rules 


2 GeV 


206 ±7 


145 ±8 


0.70 ±0.04 


[34] 


Lattice 


2 GeV 






0.72 ±0.02 


[35] 


Lattice 


2 GeV 






0.742 ±0.014 



TABLE II: AdS/QCD predictions for the vector and tensor couplings of the p meson compared to 
Sum Rules predictions, lattice predictions and experiment. We use rrif = 0.14 GeV to make these 
predictions. 



The twist-2 DAs can be expanded as [4, 5] 



-(z,fi) = 6z(l - z) 



(58) 
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where C^ 2 (£) are the Gegenbauer polynomials and £ = 2z — 1. Standard Sum Rules 
predictions are usually available only for . The twist- 2 DAs are thus approximated as 

3 



l + a!J'»-(5^-l) 



(59) 



i.e. by keeping only the first term in equation (58). We use here the Sum Rules estimates 
given in reference [36]: a| = 0.10 ± 0.05 and = 0.11 ± 0.05 . Reference [36] also gives 
explicit expressions for the twist-3 DAs: 

(3 - 30£ 2 + 35£ 4 ) . (60) 



+ 



9 II / \ 15 ||, , 15 ,|, . 
a 2W + 32 w 3l^) - ^"aW 



and 



g^\z^) = 6z(l-z) 



1 + ( \4m + + ^*) - tfZlM ) c 3 2 /2 (o 



(61) 



The Sum Rules estimates are Q(2 GeV) = 0.020 ± 0.009, w|(2 GeV) = 0.09 ± 0.03 and 
w|(2 GeV) = -0.04 ±0.02 [36]. 

In figure 2, we compare the AdS/QCD twist- 2 DAs to the SR twist- 2 DAs at a scale 
H = 2 GeV. We note that, as is the case for the AdS/QCD tensor coupling, the AdS/QCD 
DAs hardly depend on // for (j, > 1 GeV and they should be viewed as parametrizations of 
the DAs at some low scale /i ~ 1 GeV. As can be seen and as was already noted in reference 
[10], the agreement between the AdS/QCD and Sum Rules twist-2 DA for the longitudinally 
polarized meson is good. On the other hand, we note different shapes for the AdS/QCD and 
Sum Rules predictions for the twist-2 DA of the transversely polarized meson. In particular, 
we find that the AdS/QCD DA has pronounced humps near the end-points and that when 
it starts decreasing, it does so faster than the SR DA. 

In figure 3, we compare AdS/QCD twist-3 DAs to the SR twist-3 DAs at a scale /j, — 2 
GeV. The agreement between AdS/QCD and SR is quite good for both the axial vector 
DA but we note a difference between SR and AdS/QCD vector DA at the end-points: the 
AdS/QCD, unlike the SR DA, falls to zero at the end-points. 
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(a) Twist-2 DA for the longitudinally polarized p meson 





0.2 0.4 



(b) Twist-2 DA for the transversely polarized p meson 

FIG. 2: Twist-2 DAs for the p meson. Solid Red: AdS/QCD DA at fi ~ 1 GeV; Dashed Blue: 
Sum Rules DA at fi = 2 GeV. 

BRANCHING RATIOS 



We are now in a position to compute the branching ratios given by 

'II 1 \A(B° (S) ^P°J)\ 2 , (62) 



Bn(B? s) p° 7 ) 



16nMl 



M 



B 



15 




(a) Axial-vector twist-3 DA for the transversely polarized p 

meson 




z 



(b) Vector twist-3 DA for the transversely polarized p meson 

FIG. 3: Twist-3 DAs for the p meson. Solid Red: AdS/QCD DA at fi ~ 1 GeV; Dashed Blue: 
Sum Rules DA at /i = 2 GeV. 

where the amplitude A(B?, — > p°7) is given by either equation (25) for B° — > p°7 or 
equation (26) for B° — > p°7 and tb, s) is the measured lifetime of the meson [15] . 
Before presenting our predictions for the branching ratios, it is instructive to compare the 
AdS/QCD and SR predictions for the integrals I{ w2 (s p ,fi) and I| w2 (//) given by equations (8) 
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and (11) respectively. Our results are shown in table III. We note that the integrals are not 
very sensitive to the precise shape of the twist-2 DA. 



Integral 


SR 


AdS/QCD 




1.902 + 2.620i 


1.590 + 2.329z 




-6.561 + 0.030i 


-8.866 + 0.027i 


/r(o,/x) 


-6.660 


-8.989 


rtw2 
l 2 


3.330 


4.495 



TABLE III: AdS/QCD and SR predictions for the integrals If 2 and I| w2 given by equations (8) 
and (11) respectively. The SR predictions are at a scale [i = 2 GeV and the AdS/QCD prediction 
are at a scale /i ~ 1 GeV. 

We next compare the Sum Rules and the AdS/QCD predictions for the integrals i| w3 and 
f* w3 given by equations (18) and (19) respectively. Our results are shown in table IV. In 
this case, the AdS/QCD and SR predictions are drastically different. The SR DA yields 
divergent integrals for both I[ w3 and I^ 3 unlike the AdS/QCD DA which leads to finite 
results in both cases. The divergent SR integrals could be estimated by introducing an IR 
cut-off but this procedure leads to a large uncertainty in the prediction for the annihilation 
amplitude [12]. 



Integral 


SR 


AdS/QCD 




oo 


0.237(0.229) 




oo 


0.036(0.034) 



TABLE IV: Sum Rules and AdS/QCD predictions for the integrals i] w3 and I^ 3 given by equations 
(18) and (19) respectively. The SR predictions are at a scale (i = 2 GeV and the AdS/QCD 
prediction are at a scale /j, ~ 1 GeV. 

It is instructive to investigate the influence of perturbative QCD scale evolution on the 
infrared divergence encountered with the SR DA. As shown in figure 4, we evolve the SR DA 
from /j, — 1 GeV to pL — 2, 3 and 5 GeV to leading logarithmic accuracy using the evolution 
in [6]. We also show the asymptotic DA, i.e. the SR DA at /i — > oo. As can be seen, the SR 
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FIG. 4: Evolution of the twist-3 SR DA. Blue: SR DA at [i = 1 (Dot-dashed), 2 (Long-dashed), 3 
(Dot-dot-dashed) and 5 (Short-dashed) GeV; Dotted Black: Asymptotic DA; Solid Red: AdS/QCD 
DA at p ~ 1 GeV. 

DAs do not vanish at the end-points and we find that the divergence problem persists at 
scales other than 2 GeV. Also shown in figure 4 is the AdS/QCD DA which, unlike the SR 
DA, vanishes at the end-points and avoids the end-point divergences. On the other hand, 
the AdS/QCD DA lacks the perturbative evolution with the scale /i and must be viewed 
to be a parametrization of the DA at some low scale \i ~ 1 GeV. This is a shortcoming of 
the AdS/QCD DA compared to the SR DA. However, we expect the AdS/QCD DA to be 
a reasonable parametrization of the DA at the scale fi = 2 GeV relevant to the decays we 
compute here although we cannot make a strong case that it will still avoid the end-point 
divergences if its perturbative QCD evolution with the scale [i is taken into account. 

Our predictions for the branching ratio of B° — > p°7 are shown in table V. In this table, 
we show how the predictions vary with the degree of accuracy of the calculation. The pre- 
dicted branching ratio computed using the leading power amplitude at zeroth order in a s 
(i.e. equation (3)) is clearly lower than the measured value. At this level of accuracy, the 
amplitude does not depend on the DAs. The leading power amplitude becomes sensitive to 
the twist- 2 DA at first order in a s and at this level of accuracy, we find that the AdS/QCD 
and SR predictions agree with each other and with experiment. We confirm that all four 
power-suppressed annihilation contributions in B° — > p°7 are numerically small. Neverthe- 
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less, the AdS/QCD DA allows us to compute the annihilation contributions beyond leading 
power accuracy without the ambiguity due to end-point divergences encountered with the 
SR DA. At the same time, the AdS/QCD DA allows us to provide a more reliable theoretical 
estimate for the branching ratio of the decay B° — > p°7 which proceeds mainly via annihi- 
lation and cannot be reliably predicted using the SR DA due to end-point divergences [12] . 
Using the AdS/QCD DAs, we predict a branching ratio of 5.5 x 1CT 10 for this decay. This 
rare decay can be enhanced by NP [12] and it would be interesting to investigate it at the 
LHCb. 



Branching ratio (xlO 7 ) for B° p°7 



DA 


Accuracy 


SR 


AdS/QCD 


PDG 


Belle 


BaBar 


tw2 + tw3 


Lead. (a]) + Anni.[a°, (l/m b ) 2 ] 




7.67 


8.6 ±1.5 


7 o i 1.7 _i_0.9 
'• 8± 1.6 ±1.0 


q 7+2.4 ,0.6 
y - ' =t 2.2 ^o.e 


tw2 


Lead, (a*) + Anni.[a°, (l/mb)} 


7.86 


7.65 








tw2 


Leading (a]) 


7.87 


7.68 








None 


Leading (a^) 


4.76 


4.76 









TABLE V: Sum Rules and AdS/QCD predictions for the branching ratio (xlCT 7 ) of B° ->• p° 7 
using AdS/QCD or Sum Rules compared to the measurements from Belle [13], BaBar [37] and the 
average value from PDG [15]. 

CONCLUSIONS 

We have used new holographic AdS/QCD DAs for the transversely polarised p meson 
in order to compute the branching ratios for the decays B° — > p°7 and B s — > p°7 beyond 
leading power accuracy. The AdS/QCD prediction for the branching ratio of B° — > p°7 
agrees with experiment and we provide a theoretical estimate for the branching ratio of 
the rare decay B° — > p°7. We find that the AdS/QCD DAs are complementary to the 
standard SR DAs: they agree with the SR predictions to leading power accuracy and avoid 
the end-point divergence ambiguity when computing some power corrections. However, in 
its present form, the AdS/QCD DA lacks the perturbative QCD evolution and it remains 
to be seen if our conclusion remains valid if this evolution is taken into account. 
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